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O . In view of the recent diffractive dijet data from CDF run II, we critically re-evaluate the 

standard approach to the calculation of central production of dijets in quasi-elastic hadronic 

1^^ _ collisions. We find that the process is dominated by the non-perturbative region, and that even 

O ■ perturbative ingredients, such as the Sudakov form factor, are not under theoretical control. 

Comparison with data allows us to fix some of the uncertainties. Although we focus on dijets, 

our arguments apply to other high-mass central systems, such as the Higgs boson. 



^ jr. cudell@ulg.ac.be 
^alice. dechanibre@ulg.ac.be 

^oscarh@physics.nicgill.ca ; permanent address: Marianopolis College, 4873 Westmount Ave. 
Montreal, QC, Canada H3Y 1X9 
*igor. ivanov@ulg.ac.be 



Introduction 

The CDF collaboration has recently published the measurement of the cross section for exclu- 
sive dijet production \T\. These are important data, as the dijet system reaches masses Mjj 
of the order of 130 GeV, i.e. the region of mass where the Higgs boson is expected. Given 
the high centre-of-mass energy involved, ^/s ^ Mjj, the process is entirely due to pomeron 
exchange. Thus the dijets are produced by the same physical mechanism which could produce 
the Higgs boson and other rare particles. The CDF data provide a good opportunity to tune 
the calculations of quasi-elastic diffractive processes which, as we shall show, are otherwise 
plagued with severe uncertainties. 

Almost twenty years ago, Schafer, Nachtmann and Schopf [2] proposed the use of diffractive 
hadron-hadron collisions, where a high-mass central system is produced, as a way of producing 
the Higgs boson. Bjorken recognised such quasi-elastic reactions, in which the protons do not 
break despite the appearance of the high-mass system, as a "superb" channel to produce exotic 
particles ^. The first evaluation was done in the Higgs case by Bialas and Landshoff [4|. It 
relied heavily on the use of non-perturbative propagators for the gluons, which provide an 
automatic cut-off of the infrared region. That calculation was then repeated by two of us [5] , 
where we showed that by properly treating the proton form factors, one could use either 
perturbative or non-perturbative propagators. However, even with the use of form factors 
screening the long wavelengths of the exchanged gluons, one remained sensitive to the infrared 
region, as the gluons had a typical off-shellness of the order of 1 GeV. Berera and Collins [6] 
then calculated the double-pomeron jet cross sections using a perturbative QCD framework. 
They identified large absorptive corrections — corresponding to the "gap survival probability" , 
i. e. to multiple pomeron exchanges — as well as large virtual corrections — coming from the 
large difference of scales at the jet vertex, and resulting in a "Sudakov form factor" . However, 
they did not model these corrections. 

Since then the subject of diffractive dijet and Higgs production has been extensively dis- 
cussed from a variety of points of view [7]- [H]. All the calculations essentially follow the 
same pattern; the variations come in the following four ingredients: 

1. The two jets are at low rapidity, and at high transverse energy Et. They then come 
dominantly from gluon production at large transverse momentum (which we shall note 
k2) and the partonic amplitude can be calculated in perturbative QCD. Furthermore, 
the two outgoing gluons are constrained to be in a colour-singlet state, and, to prevent 
colour flow in the t channel, an extra screening gluon is exchanged (Fig. l.a). The exact 
relation between the parton k2 and the jet E^ will be discussed in Section 16.21 

2. As the gluons go from low transverse momenta in the proton to high ones in the jet 
system, there are enhanced double logarithms from vertex corrections (Fig. l.b), which 
must be resummed to give the Sudakov form factor. 

3. The infrared divergences of the gluon propagators linked to the proton are regulated via 
the use of an impact factor [5j, as in Fig. l.c, or of non-perturbative propagators [18], 
or via the Sudakov form factor linked to the hard scattering [8]. 
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Figure 1: A sketch of the various steps of the calculation. 



Finally, it is possible that the two protons interact at long distances. This leads to 
screening corrections, as sketched in Fig. l.d, and to a gap survival probability given by 
the average squared norm of the S-matrix element |(pp|>S'|pp)p jTOj [20]. 



In this paper we want to evaluate, in view of the CDF data, the uncertainties in the various 
ingredients of the calculation of central exclusive dijet production. As discussed above, ours 
is not the first such calculation, however we believe to shed light on several important points. 
We show that the calculation still lies mostly in the non-perturbative region, and that many 
standard approximations cannot be justified. In particular, we find the following: 

• At the Born level, exact transverse kinematics is important. That is to say, the momen- 
tum transfers to the hadrons cannot be neglected with respect to the momentum of the 
screening gluon, which makes the exchange colour-neutral. In Section [1], we present a 
detailed account of the part of the calculation which is under control, i.e. lowest-order 
qq — >■ qqgg via a colour-singlet exchange. There we take the exact transverse kinematics 
into account. Furthermore, diagrams in which the screening gluon participates in the 
hard sub-process cannot necessarily be neglected, as has been assumed in previous works. 
In Section 12. 3[ we show they could be important once the Sudakov suppression is taken 
into account. These issues affect the calculation of the partonic amplitudes mentioned 
above in ingredient [H 

• The leading and subleading logs that are resummed to give the Sudakov form factors 
are not dominant for the momentum range of the data, i.e. the constant terms are nu- 
merically important. This affects the calculation of double logarithm vertex corrections 
mentioned above in ingredient [2l We discuss the Sudakov form factor and the problems 
associated with the large virtual corrections in Section [2J 

• The colour neutrality of the protons has to be implemented independently of the Su- 
dakov suppression. This affects the calculations [8] that make use only of the Sudakov 
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Figure 2: Kinematic conventions for the central two-gluon production. 



form factor to regulate the infrared divergences of the gluon propagators mentioned in 
ingredient [31 This is discussed in Section [3l where we consider various ways to embed 
the perturbative calculation into a proton. 

In addition to these points, there remains the issue of gap survival, ingredient H] in the list 
above. Providing an accurate numerical estimate of this is beyond the scope of this paper but 
we do discuss the gap survival probability and current estimates in Section |H 

In view of the above uncertainties, we try to outline a few scenarios that do reproduce 
the dijet data, and that can be extended to Higgs production at the LHC After a few extra 
corrections, we give in Section [5] a simple estimate of the cross section, followed in Section [6] 
by detailed numerical results. Finally in Section [7] we present our conclusions. 

1 The lowest-order perturbative QCD calculation 

1.1 Kinematics 

The backbone of the central quasi-elastic production of two high-ii^^ jets is the partonic sub- 
process qq — > qqgg, in which the produced colour- singlet two-gluon system is separated by 
large rapidity gaps from the two scattering quarks. The kinematic conventions are shown in 
Fig. [2J We assume that the quarks are massless and consider the collision in a frame where 
the incoming quarks have no transverse momenta. Their momenta q'^ and p'^, with s = 2p ■ q, 
will be used as the lightcone vectors for the Sudakov decomposition of all other momenta en- 
tering the calculation. The momentum transfer to the first and second quarks are —ki and ^3, 
respectively, and are dominated by their transverse parts — ki and h^ (we write all transverse 
vectors in bold). The momenta of the two produced gluons are 

< = a,p^ + Ag'^ + (ki - k^)'^ , a,f3,s = {k, - k,f, 

r^ = a2p'' + (32q^ + (ks - k,Y , as/^as = (ks - ks)^ . (1) 

The largest contribution to the cross section will come from the region where longitudinal 
components obey 

l>/3i,/52>— , l>a2, «!>— , i = l, 2, 3, (2) 

s s 



as the invariant mass squared of the two-gluon system 

^aa = in + ^2) = ^-5 (3) 



,2 _ . , N2 _ [/32kl + Ak3 - ( A + /?2)k2] 



is much smaller than s. The differential cross section can then be written as a convolution 
over a phase space factorised between light-cone and transverse degrees of freedom: 

The longitudinal phase space can alternatively be written as 

d(3id(32 d(3dx r—— f3i 

— — — = -5- — , where /? = VP1P2 and x = — . (5) 

Pi P2 P X P2 

If the two gluons are both integrated in the whole available phase space, then an extra 1/2 
should be put in the expression of the cross section due to Bose statistics. 

1.2 Simplifications for the imaginary part 

As we shall see, the lowest-order calculation will lead to an amplitude which grows linearly 
with s. As the exchange is C = +1, the amplitude is then mostly imaginary, and can be 
calculated via their standard cuts. At the same lowest order, the real part is suppressed by a 
power of s, however it will be only logarithmically suppressed at higher orders. In principle, it 
can be obtained via dispersion relations, but we do not concern ourselves with its contribution, 
as it will be much smaller than the large uncertainties in the other parts of the calculation. 

In contrast to central Higgs production, two gluons can be emitted from all parts of the 
diagrams in many different ways, which is represented by the grey circle of Fig. [21 However, if 
one calculates the imaginary part of the amplitude, then there are multitudinous cancellations 
among different contributions due to the positive signature and colour-singlet nature of the 
exchange, as well as to the presence of large rapidity gaps. We show in Fig. [3] two typical 
cut diagrams that give rise to an imaginary part for rii + 77-2 jet production. The dashed line 
represents the kinematic cut of the diagram, i.e. it indicates which propagators are put on 
shell in the loop integral. The contributions of the "wrong cut" such as, for example, those 
shown in Fig. Ob cancel one another. This means that one can write the amplitude as a 
sequence of two sub-amphtudes, gauge invariant on their own: 2 ^ 2 + ni ^ (2 + 712) + ni, as 
in Fig. [3la. Since in our case ni + n2 = 2, there are three generic situations: (ni = 0, n2 = 2), 
{rii = 2, 77-2 = 0), or {rii = 1, 722 = 1). 

Each of the sub-amplitudes in Fig. Oa describes the emission of one or two gluons in the 
central region. In principle, all the vertices needed for the calculation of two-gluon produc- 
tion in arbitrary kinematics can be found in the literature. Production of one gluon can be 
described by the standard Lipatov vertex |2T], while emission of two gluons involves an effec- 
tive four-gluon vertex in the quasi- multi-Regge kinematics [22]. Such non-local vertices take 
into account gluon emission not only from the t-channel gluons themselves, but also from the 
quarks. In the case of large transverse momentum of the produced gluons, which is the focus 
of the present paper, the situation simplifies, since emission from t-channel gluons is dominant, 
and the amplitude is more conveniently calculated using Feynman diagrams and cutting rules. 
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Figure 3: Schematic representation of the imaginary part of the production amphtude of 
ni + n2 gluons in the central region, s-channel cuts of the diagrams, such as (a) contribute to 
the imaginary part, while "wrong cuts", such as (b), do not. 



1.3 Central production of two gluons with large transverse momen- 
tum 

We are interested in the quasi-elastic production of two gluons with large relative transverse 
momenta of the order of tens of GeV. The requirement that the protons remain intact effec- 
tively cuts the differential cross section at small values of momentum transfers, k^, k| < l/Bp, 
where Bp is a typical proton elastic slope in hadronic reactions (for more discussion on what 
numerical value for Bp would be most appropriate in our case, see Section [3l3!) . Therefore, 
k^, k| -C k^, and the process can be viewed as a collision of two nearly-coUinear but energetic 
gluonqj g*g* — »• gg accompanied by an additional exchange of an extra screening gluon to 
restore the neutrality of the t-channel colour exchange. The set of diagrams to be considered 
is then reduced to those of Fig. |U and to their counterparts where each gluon is emitted from 
the other side of the cut. 

Lowest-order diagrams with gluons emitted from different t-channel gluons, such as the 
one of Fig. Hlb, are suppresseqj by one extra power of Vi^- The situation may become different 
at higher orders, as we shall explain in Section 12.31 Also, as we are calculating cut diagrams, 
the sub-amplitudes are at the tree level and there are no ghost contributions. 

The Mandelstam invariants in the two-gluon collision can be written 
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A/52 



^99 



,2 I31+P2 
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(6) 



Note that these invariants depend only on the ratio x = /3i//32, which is related to the difference 
of the rapidities of two produced gluons. 

Let us also introduce the momenta of the colliding gluons 



K 



{k + k 



1^ 



KjO 



-{k + /C3 



\t^ 



(7) 



^Strictly speaking, in the lab frame, one of these two gluons can be very soft and emitted at large angle 
relative to the quark collision axis. However, after an appropriate longitudinal boost the above description 
will become true. 

^Note that the presence of hard transverse momentum in the s-channel partons in these diagrams changes 
nothing since these partons are on-shell, so that the transverse momentum does not suppress the amplitude. 




Figure 4: (a) A generic diagram for qq -^ q + gg + q scattering in the high-k2 regime; the 
shaded blob represents the diagrams for the elastic gg scattering, (b) A diagram with two 
gluons emitted from two t-channel legs is suppressed by an extra l/kg factor. 



The imaginary part of the amplitude can be represented as 

^4 6"'' r d^k 



ImM 
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where 



.(1)* 

.^Ai 



u{q + fc)7 



^Ai^l 



-K 



,■(2)* 



u(p - fc)7 ■ et m(p + fcs) 



(9) 



are the amplitudes q -^ qg with gluon polarisation vectors discussed below. Strictly speaking, 
the two colliding gluons are virtual. In our kinematics, their off-shellnesses are — /tf and — ^3, 
which are much smaller than \i^. Since the hard scale of the scattering sub-process is given 
by k^, one can neglect the non-zero virtualities and approximate the amplitude by gg — > gg 
scattering of on-shell transversely polarised gluons, and neglect the contribution from the 
longitudinal polarisations. The calculation is then manifestly gauge invariant. 

The helicity amplitudes for the tree-level scattering of two gluons in a colour-singlet state 
are most easily calculated in the centre-of-mass frame: 



M- 



AiA2^A3A4 






A3A4) 



(10) 



where (p is the azimuthal angle of the two-gluon-pro duct ion plane with respect to the quanti- 
sation axis. Note that this quantisation axis is arbitrary, and changing it will produce changes 
in M which will be compensated by opposite changes in Jx ■ 

The fact that Eq. ( ITOl) is written in a frame which is different from the laboratory frame 
does not pose any problem. Indeed, in order to pass from the laboratory frame to the gg centre- 
of-mass frame with gluons colliding along the z axis, one first has to perform a longitudinal 
boost to make the energies of the colliding gluons equal, then a transverse boost to make 
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the total momentum of two gluons zero, and then rotate the frame to ahgn the z axis with 
the direction of the incoming gluons. The large longitudinal boost does not change 0, while 
the transverse boost and the rotation by a small angle have negligible effect on the hard 
momentum k^. Therefore, one can safely understand in Eq. fITOl) as the azimuthal angle in 
the lab frame. 

The non-zero A{\\\2 — > A3A4) in Eq. ( ITOl) are 



A(++- 




J.2 


^(+-- 


-. -+) = A(-+ - 


^gg 



This list exhibits the total helicity conservation rule, which is a consequence of the helicity 
properties of a general tree-level n-gluon scattering amplitudes, see e.g. [23]. In our case it 

implies, in particular, that -|--|- and amplitudes do not interfere with any other. 

The fact that the colliding gluons are soft, and the momentum hierarchy of Eq. ([2]), sim- 
plifies the calculation. Each of the polarisation vectors for the initial gluons can be chosen 
orthogonal to both k^ and ^3, and within our accuracy can be generically written as 



2 



e^ = e^ + 



Ba • Ki ^ eA ■ K3 



/?1 + /?2 «! + 02 



:ii) 



Here, Ba is the standard polarisation vector in the transverse plane, 

eA = — ^(A, i), 

with A = Ai for the first gluon and A = — A2 for the second one, since they move in opposite 
longitudinal directions. One can now simplify 

.(1)* _ e^^-Ki _ 2 ^1 1.. |^-»Ai^i 

,lf - 2%1^ = ^^|..3|e-^^S (12) 

«! + ^2 ttl + ^2 V2 

where 0i and 03 are the azimuthal angles of k^ and ks, respectively. 
Squaring the amplitude, one finds the following expression 






.(1)*,(2)*,(1'),(2') 



Y. Y. ^C^r^, 'ri ' Ma,a.^/m:, ,, ^, , (13) 

where / labels the polarisation states of the final two-gluon system. Summation over final 
fermions and averaging over initial ones is also implied here. Note that, in contrast to the 
standard Weizsacker- Williams approximation, where the initial particle momentum is the same 
in j^ and jy, here it is different due to k 7^ k'. This induces correlations between the colliding 



gluons, and will lead in a moment to the conclusion that the fully unpolarised qq ^ qq + gg 
cross section integrated over all phase space is not proportional to the unpolarised gg -^ gg 
cross section. 

The only non-trivial interference here is between M_|_ >j and M_^ >j, with / = H — or — h. 

Such a term introduces an azimuthal dependence for high-|k2| gluons via the factor exp(4z0). 
It contributes to the azimuthal correlations between the high-i^^ jets and the proton scattering 
planes, but if integrated over (still keeping the cross section differential in kg), this term 
vanishes. This allows us to consider only the diagonal contributions in Eq. ( TT3l) . A'^ = Ai, 
A2 = A2. The result is 

f A, 

- ° \'^^\\^^\\<\K\ [|Mopcos(0i - 03 - <P[ + ^P's) + |M2pcos(0i + 03 - 0; - ^P's)] 



(/?i+/?2)'(ai + «2; 



2 



'11.^/ 



I'^iN'^sN'^iN'^a 



7/, ^fi^2(^ ^„^2 (|Mor + |M2p)cos(0i-0;)cos(03-03) (14) 

(pi + P2)n"i + «2)^ L 

+ (|Mo|' - iMsH sin(0i - 0;) sin(03 - 0;^) 

(Pl +P2)n"l + "2)^ l- 

+ (|Mop - IMsHKa^i ■ A^3)K ■ O - (a^i • 4)K • '^s)]}. (15) 

Here, |Mop and IM2P are the amplitudes squared of elastic collision of two gluons with total 
helicity A = Ai — A2 = or 2 summed over the final helicity states: 

IM2I' = i[|^+— +-l' + l^+— +l'+l^-+— +l' + l^-+-+-l']- (16) 

Note that if k were equal to k', then 0j would be equal to 0^, and one would end up with the 
unpolarised cross section, oc |Mop + IM2P in Eq. flT^ . in the spirit of the usual Weizsacker- 
Williams approximation. Another observation is that in the limit 

|ki|, |k3| < |k|, |k'|, 

the angles are 03 ~ 0i + vr and 03 ~ 0'^ + tt, so that Eq. (TT4|) would simplify to 

[|^o|' + |M2|2cos(20k-20kO] . (17) 



(A + /52)2(ai + a2)2 



After the azimuthal averaging over 0k and 0k', the IM2P term vanishes and only the amplitude 
with A = contributes, so that in this limit one obtains the Jz = rule [S]. 
Finally, since 
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Figure 5: A higher-order diagram for qq -^ Q + 99 + Q scattering leading to large doubly 
logarithmic corrections in the Feynman gauge. 



one can write the partonic differential cross section as 
da, 



1 fg^y dpidP2d^^i£i^2d^^3 f rf^k 
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where the first term in brackets corresponds to A = 0, and the second one to A = 2. 



2 The Sudakov form factor 

2.1 Resummation 

In the diagrams of Fig. |Ha, the gauge field goes from a long-distance configuration to a short- 
distance one. This is the situation for which large doubly logarithmic corrections, ~ log (k^), 
are expected, from virtual diagrams such as those of Fig. [5l These corrections are there for any 
gauge theory [21] and have been calculated in QCD in [25]. For initial gluons on-shell, these 
corrections actually diverge, and this divergence is cancelled by the bremsstrahlung diagrams 
in inclusive cross sections. Hence, the cancellation of infrared divergences means that the 
logarithmic structure of the virtual corrections is identical to that of the bremsstrahlung 

diagramCl- 

So there are two interpretations of the Sudakov form factor. On the one hand, one can 
view it as a resummation of the double-log enhanced virtual corrections. For example, at 
the one-loop order, these corrections include diagrams with the integrals over the fraction of 
the light-cone momentum, z, and the transverse momentum, q^, of the particle in the loop. 
Each of these integrals builds up a logarithm coming from regions A ^ z and £^ <^ q^ ^ /i^ , 
respectively. Here, A is a cut-off to be discussed below, P is a typical virtuality of the initial 



^ For the case at hand, the infrared region is cut off by the ofF-shellness of the initial gluon, so that the 
logarithms are large, finite and of the order of log(sgg/(ki +k)^) log(sgg/(k3 +k)^) > 25 for the CDF exclusive 
jet production. 



gluons in the g*g* —>■ gg process, and fi is the scale of the hard process, which is of the order 
of the jet transverse energy Et- 

On the other hand, the cancellation of logarithms between virtual and real diagrams means 
that same Sudakov form factor can be reworded in a Monte- Carlo language as the probabil- 
ity of not emitting any extra partons. In this interpretation, one considers collision of two 
gluons and calculates the probabilities that a hard sub-process is accompanied by a certain 
number of secondary partons. Subtracting from unity all the probabilities of emission gives 
the probability of the purely exclusive reaction. 

Using this technique, one can then use for the Sudakov form factor the expression coming 
from Monte Carlo simulations [26] : 

T(/i^ f) = exp [-s(/i^ f)] , ^(/i^ f) = f ^^fls!) f ^ dz [zp,, + NfP,,] . (19) 

j(^2 q MX Jq 

Here, the lower scale i'^ is understood as the virtuality from which the evolution starts, and 
Pgg and Pgq are the unregularised splitting functions. 



Pgg{z) = 2N, 



z 1 — z .^ , 
+ + z(l -z) 



Pg,{z) = \[z' + {l- zf] . (20) 



We shall refer to this form of the vertex corrections as the Splitting Function Approxima- 
tion (SFA). 

If instead of as(q^), we take as(At^) in ( p^Qj) . we can easily work out the double logarithm 
approximation (DLA), which comes from the 2A^c/(l — z) term in Pgg-. 



2 „2^ 3a, /-^ f^q ,_„/ 1 



s,u,',n = ^]^^io,y^). (21) 

The coefficient in front of the double logarithm depends on the definition of A. First, pure 
DGLAP kinematics leads to the cut-off A = q^//i^. Using this cut-off, one obtains the expres- 
sion worked out in |25j. However, it can be argued that the coherent effects lead to angular 
ordering in the successive splitting of the secondary partons [27]. This ordering introduces a 
more restrictive limit on the z integral, which becomes linear in |q|: A = |q|//i. With this 
definition of A, one obtains a double-log result, which is twice smaller than in [25] : 

S...(,.^<''^) = |;^log^(^). (22) 

Numerically, this change is very substantial; it can easily lead to factors (9(10) in the cross 
section. In the following, we shall only use the latter prescription for A. Note that, together 
with formula flTOl) . it can also be derived from the CCFM equation [28] . 

Several comments are in order. First of all, Eq. flT9l) does not include all the single loga- 
rithms present in the process, which, to the best of our knowledge, have been evaluated only 
for the quark vertex |29j. One does not know in the present case whether further single- 
logarithmic corrections are large, or whether they break the exponentiation. Moreover, even 
if one takes Eq. (TT9l) at face value, its validity depends on the fact that the logarithms are 
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Figure 6: The contribution of the various terms entering the argument of the Sudakov form 
factor. 



dominant. We show in Fig. [HI the contribution to S of the constant terms, compared to that of 
logarithms, for /i = 10 GeV. For typical scales relevant for the CDF measurements {fi ~ (9(10 
GeV), i ^ 1 GeV) the single-log corrections are large. They are negative for Pgg and positive 
for Pqg. One can even say that the overall contribution oi g ^ qq splitting is equally important 
as g —>■ gg. As we explained above, these logarithmic corrections are only an educated guess, 
and the fact that they are huge makes the theoretical predictions unstable. 

Furthermore, as can also be seen from Fig. [6l it turns out that even the constant terms in 
flTOj) are numerically important. One sees that, for an upper scale of 10 GeV, the logarithms are 
dominant only for i -C 2.5 GeV, which means that we can trust the perturbative formula ( I19p 
only in the non-perturbative region! 

Our analysis implies that any other source of single-log contributions (for example, the one 
discussed in Section IX^ below, or a mere shift of the upper and lower scales in the logarithm) is 
also expected to affect strongly the resulting value of the Sudakov form factor. In Section [6l3] 
we will study this sensitivity in detail. 



2.2 What is the upper scale in the Sudakov integral? 

The standard discussion of the Sudakov form factors relies almost exclusively on corrections 
to the quark electromagnetic form factor. In this case, the hard vertex is characterised by the 
single kinematical quantity Sqq. In the Monte Carlo language, this energy defines the phase 
space available for the secondary parton emission which needs to be suppressed. This is why 
the upper scale in the Sudakov integral can be taken as /i^ = c ■ Sqq, with numerical coefficient 
c = 0{l). 

The same argument applies to the Higgs central exclusive production. In this case, the 
hard vertex is effectively point-like, as the transverse momenta and virtualities of the top 
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quarks inside the loop are larger or of the order of the Higgs mass. 

With dijets with large invariant mass, we enter a new kinematical regime. The two gluons 
can have large invariant mass Sgg both via large transverse momentum or via strong longitudi- 
nal ordering, see ([3]). In the former case, the transverse momentum exchange in the gg —>■ gg 
scattering is large, k^ ~ Sgg, while in the latter case, kg <C Sgg. 

A typical double-log-enhanced correction to some vertex V gets one logarithm from the 
fraction of the lightcone variable and another one from the transverse momentum or virtuality 
integral. In order for the transverse momentum integral to produce a sizable log, there should 
exist a large transverse momentum or large virtuality inside the effective vertex. In this case 
the transverse momentum in the loop, when flowing through V does not change its amplitude. 
In the Monte Carlo language, in order for the backwards evolution to develop, one must 
have a large initial hard virtuality inside the vertex. On the contrary, if the vertex does not 
involve any large transverse momentum or large virtuality, then the transverse loop integral 
is suppressed by the vertex V, and the transverse logarithm does not build up. Once the 
structure of the vertex V is resolved by the incoming gluons, the structure of the amplitude 
changes, and double logarithms disappear. Because this transition is not sharp, the actual 
value of the upper scale is somewhat uncertain, but it must be of the order of k^. 

This can be immediately seen in the extreme case of production of two gluons in multi- 
Regge kinematics (and therefore large Sgg) but with small transverse momenta, of the order 
of typical loop transverse momentum in the BFKL ladder. The amplitude of the gg -^ gg 
subprocess is then oc Sgg/tgg, and this process does not involve any hard gluon. Consider 
now a loop correction to it. If the loop integral involves large transverse momentum q, then 
it will flow through the "original" gg —>■ gg subprocess and will suppress its amplitude to 
oc Sgg/c^. This suppression prevents the development of the transverse momentum logarithm, 
in accordance with the general BFKL theory. It is in this suppression of the vertex where the 
key difference lies between the usual Sudakov correction to the quark electromagnetic form 
factor with its point-like vertex (or the Higgs production) and the dijet production. 

With this discussion in mind, we believe that the physically motivated choice for the upper 
scale in the Sudakov integral should be related with the transverse momentum transfer in the 
99 ~^ 99 subprocess, but not with Sgg. 

As for the present study, with dijets produced in the central rapidity region and not in 
multi-Regge kinematics, these two prescriptions do not lead to any major difference in the 
Et, rain distribution of the cross section. However, they will lead to vastly different M'jj shapes 
of the cross section, which will be discussed in Section 16.21 

2.3 The role of the screening gluon 

It is usually assumed that the extra gluon in the t-channel is needed only to screen the colour 
exchange and does not participate in the hard sub-process. Fig. HI a. So the pomeron fusion 
is very similar to gluon fusion (up to colour factors), and the details of the final state of the 
hard interaction do not substantially change the calculation. This would make dijet exclusive 
production essentially identical to Higgs exclusive production, and leads to the hope that the 
CDF data can be used to calibrate Higgs production. At the lowest order, this assumption is 
justified, since the emission of two gluons from two different t-channel legs. Fig. Hlb, leads to 
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Figure 7: At the Born level, (a), the diagram with two gluons emitted from different legs 
involves hard momentum flow in the t-channel, represented by thick gluon lines. As a result, 
the one-loop corrections, (b), are never double- log enhanced. However, they "localise" the 
hard momentum flow in the diagram, and successive loop corrections are now double-log 
enhanced. Their calculation amounts to calculation of a novel type of Sudakov form factor to 
the irreducible gggg -^ gg process, (c). 



an extra l/k2 suppression. 

The virtual corrections discussed so far in connection with the Sudakov form factor do not 
lead to any non-trivial "cross-talk" between the two t-channel gluons. Besides, the standard 
BFKL-type exchanges can be re-absorbed in the evolution of the gluon density, as we shall see 
in the next section, so that the screening gluon effectively decouples from the hard sub-process. 

Here, we would like to discuss two potential mechanisms through which the screening gluon 
may get involved in the dynamics. 

The first one is specific to dijet production. We know that the Born-level amplitude of 
the standard diagram Fig. HI a must be corrected by a Sudakov form factor. Keeping only the 
leading powers of the transverse momenta, we can write 



ImM ~ 1^/ (^k^exp [SiklX)] ~ ^^^exp [-S{kl (k^))] 



(23) 



where (k^) ^ (9(1 GeV^) is the position of the saddle point. 

Let us now consider similar corrections to the suppressed diagram of Fig. |Hb, the hard 
part of which is shown in Fig. Ha. At the one-loop level, corrections to these diagrams are 
not double-log enhanceqj. One can however channel the hard transverse momentum flow as 
shown in Fig. [71b. This as correction is enhanced by a single log. After that, however, the 
lines corresponding to the hard process remain the same at higher orders, and the outer gluon 
legs do not carry any hard momentum, see Fig. [71c. Thus, successive loops do produce double- 
log corrections to the irreducible gggg -^ gg vertex, and, perhaps, may even be resummed, 
leading to a novel type of Sudakov form factor, the calculation of which might be interesting 



*We remind the reader of the QED result for the Sudakov form factor in qq coUision with quark virtuahties 
ki and fc| and total momentum q, [24ll30j : 



27r 



log 



log 



In a very asymmetric case, with pi ^ q^ and P2 ^ q^, this expression is only single-log enhanced. 
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(a) 



(b) 



(c) 



Figure 8: A colourless final state / can be produced not only via the standard diagram (a) 
but also via a collision of two gluons produced in multi-Regge-kinematics, (b) and (c). The 
diagram (b) is the first term in the Sudakov form factor, while diagram (c) is specific for 
central exclusive production and modifies this form factor in the subleading logarithms. Note 
that the hard process is totally concentrated in the blob and that no gluon line carries large 
transverse momentum. 



on its own. However, we stress that the corrections to the lowest-order result always lack one 
logarithm, as they are of order as log(k2) {as log^(k2))". 

Arguably, this is an indication that the corresponding factor that accompanies the Born- 
level amplitude, which we write as exp[— S'„e«,(k2, k^)], is not as small as the Sudakov form 
factor: it might be that S'„e^(k2,k^) -C S'(k2,k^). A very rough estimate of the resulting 
amplitude is (see derivation in Appendix A) 



ImMr 



-^ J -^ ^^P [-Sneu,{\^iy)\ . 



(24) 



Note that in contrast to the "standard" amplitude (l23l) . the integrand here extends to much 
higher values of k^. Therefore, it might happen that, after all, the diagram Fig. HJb is not 
as much suppressed as it looks at the Born-level. Certainly, for a very hard process, i.e. for 
k2 — > cxD, it can be safely neglected. However, its importance grows at smaller E^, and it is 
not obvious without a detailed calculation from what values of E^ the estimates from Eq. (1251) 
and Eq. fl2^ become of the same order, and whether this interval includes the CDF kinematic 
region. As the t-channel gluons are in a colour-singlet state, the diagrams of Figs. H^ and b 
have opposite signs, so that the overall effect will be to decrease the jet cross section w.r.t. 
the Higgs or 77 production cross sections. 

The above corrections are specific to the dijets and are absent for Higgs or 77 production. 
Let us now discuss another potential correction that affects all of these final states (which we 
denote generically /) and, even more importantly, which is not suppressed by a power of the 
hard scale as was fl24|) . 

If a process gg —^ f is possible, then one can produce / not only via the standard mech- 
anism. Fig. [Hla, but also via processes Fig. [Sib and Fig. [Sic. Note that none of the gluons 
shown in these diagrams carries hard transverse momentum of order Ej-. Diagrams Fig. [Hlb 
and Fig. [Hlc possess an extra as but are enhanced by at least one large logarithm. More specif- 
ically, diagram Fig. [Hlb is double-log enhanced and is effectively taken into account by the 
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Sudakov form factor. However, the standard Sudakov form factor (the one used in inclusive 
production of a colourless state) does not include diagrams like Fig. [Hlc, which are specific for 
the central exclusive production. 

We stress that the contribution of this diagram is not small, as it is single-log enhanced (see 
Appendix B). Therefore, such diagrams must be resummed, and they will lead to modifications 
to the Sudakov form factor at the single-log level, that were absent in the inclusive case. 
Recalling the very substantial role of single-log effects in the Sudakov form factor, one might 
expect numerically very sizable corrections due to the diagrams like Fig. [Hlc. Further analysis 
is definitely needed to bring these corrections under control. 

3 Embedding the gluon production into proton-proton 
collisions 

3.1 Impact factors 

So far, we have considered colour-singlet quark-quark scattering, and expression (ITSll is singular 
when the exchanged gluons go on-shell, i.e. k^K^Kg 'kl'^n'^n'^ — * 0. It has been argued [8] that 
this divergence could be regularised via the Sudakov form factor, provided one chooses the 
smallest gluon momentum as a lower bound for the q^ integral in (IT9l) . This is in general 
not sufficient, as the same diagrams lead to a contribution to the inclusive jet production 
cross section, which will not contain any Sudakov form factor, but must be finite nevertheless. 
Furthermore, is well-known [31] that the Sudakov form factor is not required to regularise the 
divergence: indeed, gluons with a long wavelength A average the colour of the proton, and 
hence a suppression proportional to Tc/A must always be present when coupling to a colour 
singlet of size r^.- This is usually taken into account via the introduction of an impact factor 
$(ka, kb) such that 

$(0,kfe)=$(k„,0) = 0. (25) 

Here k^ and k^ are the transverse momenta of the two t-channel gluons that couple to the 
proton. The differential cross section for two-gluon production in proton-proton scattering 
can then be represented as 



d(Jpp = d(Tqq (g) J]^ 



=1,3 I- 



T(ki,/^2)$(k,-.v, 



T(ki,<2)$(k',-^0 



(26) 



The symbol ® here indicates that the factors are to be introduced inside the loop momentum 
integrals / (P'k cP'k' . 

The presence of the impact factors removes all the infrared singularities present in the 
partonic- level expressions. Let us also note that g^"^ in ( TT8l) should be understood as gfia.rd dsoft^ 
where ghard is a hard coupling coming from gluon-gluon collision, while gsoft is the coupling 
at the gluon-proton vertex. 

Note that when linking hadronic and partonic cross sections as in (|26ll . we do not need 
to worry about the flux factors and the phase space. In principle, one could write the cross 
section for proton- antiproton scattering from the start, which would be expressed in terms of 
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the hadronic fluxes and phase space, and then work out the two-gluon production amphtude 
in the proton-antiproton colhsion. The net resuh would consist in the replacement of s by Spp 
both in the flux/phase space and in the amplitude, and since the lowest order cross section 
is independent of the energy, these changes would cancel each other. This is consistent with 
the fact that the partonic lowest-order cross section flTSj) depends on the momentum fractions 
Pi and /?2 only via there ratio x = /3i//32, so that j3i can be taken either with respect to the 
quark or to the (anti)proton. 

This situation changes beyond the lowest order. In particular, we shall introduce below the 
unintegrated gluon density of the (anti)proton that depends on the fraction of the hadronic 
momentum carried by the gluon. Therefore, from now on, we shall understand the fractions 
Pi and Oj with respect to the hadrons, i.e. as if the decomposition (Q involved the hadronic 
momenta P^ and Q^, rather than the partonic momenta p^ and g'^. 

3.2 Form factors from a quark model 

The simplest approach to the quark form factors is to consider that the ultra-relativistic proton 
is dominated by a 3-quark Fock state [32l [33] . One can then derive the form factors in terms 
of the quark light-cone wave function, and show that the form factor corresponding to the 
two gluons of momenta k^, kf, coupling to the same quark is identical to the measured Dirac 
helicity-conserving form factor Fi ((k^ -|- k^)^). The contribution of the couplings to different 
quarks can then be parametrised as Fi (k^ + k^ + cka ■ kf,): 

<l>(k„ k,) = 3 [FiiiK + kbf) - Fi (k^ + \4 + cK ■ kb)] , (27) 

This makes the jet production finite. An s dependence can be easily introduced, provided one 
assumes that a simple-pole pomeron, of intercept Oq and slope a', is a good approximation. 
One then has to multiply the form factor $ by the Regge factor [3] 

Rjets = l/C"""'* (28) 

with t the total momentum carried by the two gluons, and C,i the longitudinal momentum 
losses of the proton or of the antiproton, in our case /3i -|- P2 or ai + 0^2. 

The coefficient c can then be fixed, as well as the soft coupling gsoft, so that one reproduces 
soft data such as the total cross section and the total elastic cross section, for which the Regge 
factor has to be changed to 

ao-\-a. t 

Rsoft = a's^. (29) 

For the CDF cuts, one has ^i ~ 0.01, so that the Regge factor of the jet cross section is 
comparable to that of elastic scattering at 50 GeV. The t-slope of the exchange is thus of the 
order of 4 to 5 GeV"^. 

The two main disadvantages of this method come from the fact that, even if simple-pole 
exchanges dominate the soft amplitudes up to the Tevatron energy, they will receive non- 
negligible corrections at LHC energies. Furthermore, the use of perturbative 2-gluon exchange 
to reproduce elastic cross sections does not work, in the sense that it produces a very large 
curvature for da/dt. So the normalisation of c and Qsoft in this approach is at best tentative. 

It may thus make sense to use DIS data to normalise the impact factors. This is precisely 
the idea behind the use of unintegrated gluon densities in this calculation. 
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3.3 Unintegrated gluon density 

The distribution of gluons inside an ultra-relativistic proton can be modelled in a diagonal 
process by replacing the Born-level impact factor of the quark model with an unintegrated 
gluon density [34| : 

CF^^k, -k) = J^BoUxg, k) ^ J^{xg, k) . (30) 

In numerical calculations we used parametrisations of the unintegrated gluon density developed 
in [3H [35] . These parametrisations were obtained by fitting, in a fc^-factorisation approach, the 



proton structure function F2p to the HERA data in the region of photon virtuality < Q^ < 35 
GeV^ and for Bjorken-x smaller than 10^^. Although the fct-factorisation approach is designed 
mostly to work in the region of small x and moderate Q^, it was verified that the inclusion of a 
valence quark contribution extended these fits of F2p data to a much broader kinematic region, 
for Q^ up to 800 GeV^ and for x up to ~ 0.5. These fits, without any readjustment, were also 
found to give a good description of the charm contribution to F2p, Fg^, of the longitudinal 
structure function, F^, as well as of diffractive vector meson production [35l [36] . This serves 
as an important cross-check of the universality of the unintegrated distributions. 

These fits of T{xg^ k) were constructed as sums of two terms, a hard part and a soft 
part, with a smooth interpolation between them. The hard component describes the effects of 
hard perturbative gluon exchange, and therefore is based on direct differentiation of conven- 
tional gluon densities (for which LO fits by GRV^S], MRS [39], and CTEQ[40j were used) and 
smoothing out the saw-like behaviour of the result. The soft part describes soft colour-singlet 
exchange in the non-perturbative regime, and is constructed in a phenomenological way in- 
spired by the dipole form factor. We stress that in the soft regime the term "unintegrated 
gluon density" must be understood simply as the Fourier transform of the colour dipole cross 
section. 

One also should not forget that at small x there is no strong hierarchy among successive 
gluons in a t-channel gluon ladder, which implies that the boundary between soft and hard 
interactions becomes very smooth. For example, because of this soft-to-hard diffusion, the 
robust feature of all the parametrisations obtained in [311 [35] was that at a; ~ 10~^, the struc- 
ture function ¥2^ received dominant contributions from the soft part for photon virtualities 
up to Q^ ~ 10 GeV^. Since the process we consider in this paper also takes place at small 
Xg, this brings further concerns on the validity of the often-cited statement that the central 
exclusive diffractive production is in the perturbative regime and is well under control. 

One of the manifestations of the soft dominance in this kinematic region is the observation 
that exponent A which controls the energy growth of the integrated gluon density (often called 
the effective Pomeron intercept) calculated within these fits is significantly smaller than the 
one calculated directly from the conventional gluon densities. At typical gluon transverse 
momenta of 1-2 GeV^, we find A ~ 0.1 — 0.2, instead of 0.3 — 0.4, as obtained from the 
conventional gluon densities. This is expected since the standard DGLAP approximation 
does not take into account new portions of the phase space that open up at small x, and 
instead attributes an artificially fast growth rate to the gluon density itself. 

The structure function ¥2^ allows one to obtain the fits of the forward unintegrated 
gluon density, while the process we study here contains skewed unintegrated gluon densi- 
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ties at non-zero momentum transfer J-'{xi, X2, k^, k;,). We construct the latter in the following 
way. The effect of skewness is effectively taken into account by assuming that, for xi <^ X2, 
jF(xi, X2, ka, kf,) behaves as the forward density taken at Xg = 0.41x2. In our case, X2 = Pi +(32 
for the upper and X2 = ai +0:2 for the lower proton. The coefficient 0.41 effectively takes into 
account the skewness factor, introduced in [!41j and used in [8J. Numerically, this correction 
amounts to a factor 0.41"^'^ ^ 1.2^ ^ 2. Calculation of this factor using the conventional 
gluon densities would give a much larger value, 1.4^ ^ 4. 

The non-zero transverse momentum transfer was introduced via a universal exponential 
factor constructed in such a way that it respects conditions ( !25|) and takes shrinkage into 
account, similarly to [35l [36] . For example, the full expression for the upper proton used in 
our calculations reads: 



1 9 9 

^(Xi,X2,k,A^l) = ^(0.41(/3i+/?2),^^ 



2k2ft? r 1 



^°+^"'^°^(at^)^^ 



(31) 



with Bo = 4: GeV-2, a' = 0.25 GeV-^ and Xq = 3.4 ■ 10"^ With these values of the diffractive 
cone parameters, the slope of the k^ distribution da/dk.1 oc exp(— i?pk^) is approximately 
equal to i?p ~ 4 GeV^^. This is significantly smaller than the slope in the elastic pp collision 
at the Tevatron energy, 

dapp{el.) 



oc exp(-25ei|t|) , Bel ~ 8.5 GeV 



-2 



d\t\ 

In our opinion, the choice Bp ^ 4 GeV~^ is more natural than Bp ^ Bf,i since the elastic 
scattering at the Tevatron energy probes gluon distribution at x ~ 10"^, while in our process 
the energetic gluons carry ~ 10~^ of the proton momenta. However we stress that this choice 
is model dependent, and it introduces an extra uncertainty into theoretical calculation of the 
central production cross sections. 

Finally, note that we incorporate the unintegrated gluon density and the Sudakov form fac- 
tor as independent factors. It might be argued [28l US] that a more correct procedure would 
be to define the unintegrated gluon density via a derivative of conventional gluon density 
times the square root of the Sudakov form factor. Either way one obtains only a convenient 
parametrisation of the true unintegrated gluon distribution function. The only essential re- 
quirement is that the same prescription be used for all the processes. We believe that the way 
our fits were constructed in [31] and are implemented here, this requirement is satisfied. 

4 Gap survival 

There is one final aspect that has to be tackled to finish the calculation, and it has to do with 
gap survival probability. This concept was introduced by Bjorken a long time ago [3J, and has 
been recently investigated in detail in [20j. We remind the reader of the main argument. As 
shown in Fig. [1], the process that we have calculated at the proton level may still have to be 
corrected for initial and final-state interactions. The argument is that, due to the fact that 
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the hard interaction occurs at short distance, and does not change the quantum numbers of 
the protons, it does not influence the rescatterings. Of course, these can change the transverse 
momenta of the protons, and one would have to convolute the hard scattering with the multiple 
exchanges. However, these correlations disappear once one works in impact parameter space, 
so that the total interaction probability can be thought of as the product of the hard scattering 
probability multiplied by the probability for the two protons to go through each other, i.e. 
the S matrix element squared |S'(b)p = \{pp\S\pp)\'^ . It is easy to work out the square of 
the absolute value of this correction from the expressions of the total and of the elastic cross 
sections. 

Starting with 

^ = T7r^K^'^)l' (32) 

at lOTTS^ 

one can use the usual definition of 5* 

5(b) = l + m(s,b) (33) 

to get the partial wave 

f (fA a{s,t) ,„^, 

a(s,h) = / - — -t: (34) 

which leads to the expressions 

atot = 2 fd^hlma{s,h), (35) 

Cei = / c^^b |a(s, b)| . (36) 

The square of the S-matrix density is then the square of the deviation of a{s, b) from i [191 EO] : 

\S{h)\^ = \t-a\^ (37) 

Hence any fit of the differential elastic cross section can be used to estimate the gap survival 
probability. Generically, the gap survival will tend to 1 at large b and be suppressed at small 
b. All present models agree on the fact that, at the Tevatron, the elastic amplitude approaches 
the black-disk limit a{s, h) = i for small b. 

The simplest fit of the elastic cross section comes from CDF, who fit their data to ^3] 

^ = iV exp {2B,it) (38) 

with N = 334.6 ± 18.8 mb GeV-^ and 2Bei = 16.98 ± 0.24 GeV^^. From this, one can get an 
estimate of the gap survival probability, by assuming that the amplitude is purely imaginary. 
One then gets 

"(^•"'"^^'^K^i^) ='(°-974±°M2)exp(-^) (39) 
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Clearly, this is very close to the black-disk limit at 6 = 0. So one may expect a substantial 
suppression of the cross section due to screening corrections. If the cross section is at really 
short distance, then the average gap survival would be around 15(0)1^, i.e. less than 0.5 %. 

Fortunately, this estimate is overly pessimistic as the cross section is not concentrated at 
very short distance. Moreover, our estimate fl39l) gives only a lower bound, the amplitude has 
a real part, and also as the measurements of the E811 collaboration |44] suggest one may be 
further from the black disk limit. Furthermore, although the differential cross section is close 
to an exponential near t = 0, it has another structure at higher t. Putting all these ingredients 
together is beyond the scope of this paper. Let us simply mention that several estimates are 
present in the literature [15], and that most of them range from 5% to 15% at the Tevatron, 
and about a factor 2 lower at the LHC 

We want however to point out a few problems with the standard calculations of gap 
survival: 

• the conjugate variable of b is ki + ks, the relative momentum of p and p. This is a 
combination of soft momenta, and folding in the gap survival (which is small at small b) 
will further shift these momenta to the long-distance region. It is then very unlikely that 
the screening corrections will be given by the simple gap survival formalism. In fact the 
screening corrections will probably have a smaller effect, as they do not suppress truly 
soft cross sections very much at the Tevatron. 

• most estimates [45j assume that the ki and ka dependences factorise. This makes the 
calculation of the gap survival quite simple. However, because of the k and k' integra- 
tions, this is not true in our case. Given the large uncertainties in the gap survival, this 
may not be a crucial issue. 

5 Rough Estimate of the Cross Section 

5.1 Bare cross section 

Before presenting detailed numerical results, we find it useful to make some simple order-of- 
magnitude estimates. 

Let us start by estimating the bare cross section at the hadronic level, i.e. by keeping 
the proton form factors but omitting the Sudakov form factor altogether. Consider again the 
cross sections ( ITSl) . (!26|) . Since the strongest dependence on the ki and ks comes from the 
proton form factors, it is reasonable to take Ki ?^ ^3 ?« k and n'l k. n'^K.'k! in the numerator. 
Assuming the correlation between k and k' to be weak, one can neglect the total helicity 2 
amplitudes, so that the cross section simplifies to 

with ahard = 9 hard/ ^''^ ^^^ Correspondingly for asoft- We introduced here a dimensionless 
quantity 



{<!>^)^l-^ld'k,d'k. 



d'^k $(k, -k - ki)$(-k, k + kg 
(k + ki)2(k + k3)2 
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(41) 



Note that the x-integration in (l40l) spans over the x > 1 region, which allows us not to include 
the 1/2 factor due to the Bose statistics of gluons. The longitudinal phase space integration 
gives 

L = log(/5max//5min) log X^ax • (42) 

The cross section integrated over k^ > k^ „j„ is 



^ ~ P ^ ^ «Lrd«so/t ('^^)> 



2,min 



97r 



which gives for |k2^. 



10 GeV 



a^O.UfihLal^.^alft i^^)- 



(43) 



(44) 



Inserting the cuts used by Tevatron, we get L ^ C(5), while a^^rd ~ 0.04. Then, one has 
to estimate ($^), which is an intrinsically soft quantity. Let us recall the expression of the 
elastic pp scattering in the same approximation: 



0"ei 



4«L/t 
8l7r2 



rf^kl 



rf2k$2(k,-k-ki) 
k2(k + ki)2 



(45) 



Assuming that the strongest ki dependence comes from the exponential diffractive factor in 
$ and neglecting the energy dependence of $, one can roughly estimate 
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(46) 



(47) 



The estimate (1431) now reads: 
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For |k2|mm = 10 GeV, it gives very roughly a ~ (9(4 fib). 

This estimate, which normalises the jet cross section to the elastic cross section, includes 
implicitly a gap survival probability, which we assume here to be of the same order of magni- 
tude for both processes. 

5.2 Sudakov suppression 

CDF has measured the dijet central exclusive cross section to be about 1 nb at E^^^ = 10 GeV, 
which is three orders of magnitude below the above estimate. It indicates that the Sudakov 
suppression indeed plays a crucial role in this process. 

The Sudakov form factor enters the loop integral, and it reshapes the k^ regions that 
contribute most to the amplitude. Before the introduction of the Sudakov form factor, the 
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loop was dominated by the soft momenta due to 1/k^ factor. Now, the weight of the soft 
momenta is suppressed in Eq. (H6ll : 



As a result, the dominant k^-region shifts towards harder scales. Using the saddle point 
approximation, it has been estimated [8] that the dominant region is at |k| ^ 1-2 GeV. It is 
often claimed that this makes the loop sufficiently hard to justify the applicability of pQCD 
and the usage of perturbative fits to the gluon density. 

Here, we discuss this issue in some detail. With two simple estimates, we will show below 
that the overall suppression and the |k| shift depend strongly on the details of the Sudakov 
form factor. 

Let us switch to a more convenient variable x = log(k^/A^), where A = Aqcd and rewrite 
(|49]) as 

-^ f dx e-^(^) , X{x) = a; - 2 log $(x) + S{x2, x) , (50) 

where X2 = log (kg /A^). 

In the first approximation, the proton form factor $ plays the role of the infrared cut-off of 
the above integral: x > 0. Therefore, without the Sudakov form factor altogether, the integral 
in (1501) is of order one. 

Now, let us take the Sudakov form factor in the double-log approximation fl22]) with the 
fixed ^^(kg). Then, 

I{x)=x + - — (a;2-a;)^ (51) 

P0X2 

and the position of the saddle point Xq is at Xq = ^2(1 — Po/Q) < 0, i.e. the integral is still 
saturated in the soft region. The estimate of the suppression factor is then 



/ rfxexp[-J(0)- J'(0)x] 
Jo 



exp [-3x2//3o 
1-^ 

^ /3o 



which is about 0.2 for k2 = 10 GeV. If, instead, we take a running ^^(k^), then in the same 
double-log approximation 

Z[x) = X -F -- ( X2 log X2 + X ) , (52) 

Po V X / 

which gives now 

■^ "•" 6 

which now makes typical k^ ~ 1 GeV. The integral then becomes 

exp — 3- log ( i + — ) , (54) 



7r/5oX2 



3(l + f)2 



/?o V 6 
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which is about 8 ■ 10^^ for k2 = 10 GeV. 

Such a severe quahtative change of the resuh clearly indicates its sensitivity to different 
prescriptions used in the Sudakov form factor. In the view of additional single-log correc- 
tions coming from unspecified scales in the logarithms as well as from the extra diagrams 
discussed in Section I2.3[ which are not yet under control, one must conclude that the claims 
of perturbativity of the k-loop are unjustified. Depending on the assumptions used, the Su- 
dakov suppression can differ by one order of magnitude, and mean k^ can vary from (9(0.1) 
to 0{1) GeV^. This uncertainty affects not only numerical results, but also some qualitative 
arguments. 

Nevertheless, we see that it is possible to get a suppression factor of the order of 100 from 
the virtual corrections. This reduces our estimate to about 10 nb. An extra suppression of a 
factor 3 comes from the shift in Et when one goes from partons to jets. We shall discuss this 
in the next section. So we see that all the ingredients of this calculation can lead to a cross 
section in rough agreement with the data. In fact, we shall show that it is possible to get 
an exact agreement with the data, but that the number of adjustable theoretical corrections 
allows for many theoretical possibilities. Hence the CDF data turn out to be very important 
for tuning the theory. 

6 Numerical results 
6.1 Cuts 

The two gluons produced through pomeron exchange hadronise into jets. As is unavoidable, 
this hadronisation can result in any number of jets, not just two. Theoretically, the difference 
between the two gluon cross section and the dijet cross section {i.e. the 3-jet veto) is a 
correction of the order of asiE^) [36]. Given the large theoretical uncertainties in other parts 
of the calculation, we do not take this correction into account. 

On the other hand, the transition from the partonic to the hadronic level involves correc- 
tions which are larger, and which are due to radiation outside the jet-finding cone — we shall 
generically refer to them as splash-out. The structure of these corrections is known [171 HH] 
and involves a constant shift in Et of the order of 1 GeV, due to hadronisation, as well as 
a correction proportional to as{ET)ET due to radiation. The corresponding shift in Et has 
been estimated, for the cone algorithm used by CDF [39] to be 



E^:;' = (0.75 - 0.80)k^^''"^ (55) 

We shall also consider a previous parametrisation [50], where 



pjet _ iparton 

J-L/rp rurp 



1 - ^ «5 ((fcr*° 



1 GeV. (56) 



One can then compare our results to the data, for which the cuts are summarised in Table [TJ 

We show in Fig. [9lb the effect of these corrections. We see that for the Et range of the 

CDF data, the effect of splash-out is non-negligiblqj and amounts to a correction of the order 



^We thank V. Khoze for pointing this out to us. 
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Figure 9: (a) Our reference curve, corresponding to the parameters of Table [21 chosen so that 
it goes through the CDF run II data [T] . (b) The suppression of the cross section for various 
splash-outs discussed in the text. The plain curve is as in (a). 



ai + 012 


0.03, 0.08 


.,« _ 1 log./' ft") i-\ 1 


[-2.5,2.5] 


Vp-ytetl^ i = 1, 2 


>3.6 


yp-yfetl^ i = 1, 2 


>3.6 



Table 1: Experimental cuts for CDF run II data [Ij. 



24 



parameter 


value 


equation 




0.20 GeV 




scale of as in partonic cross section 


■^99 




scale of as in Sudakov form factor 


q^ 


m 


angular ordering 


yes 


A = |q //i in ([19]) 


terms in Sudakov exponentiation 


constant+log+log^ 


m 


lower scale of Sudakov integral 


£2 = (k + k,)2 


m 


upper scale of Sudakov integral 


/x^ = ki/2 


m 


unintegrated structure function 


ref. [36J 


m 


gap survival probability 


(^2) = 15% 


is]) 


splash-out 


Ef * = 0.8£;^'^^*°"^ 


m 



Table 2: Default parameters of the reference curve of Fig. [9l 



of a factor of 3 for EJp"^"^ = 10 GeV, and the various possibilities for the splash-out bring in 
an uncertainty of the order of 1.7 at E^^" = 10 GeV to 4 at E^^" = 35 GeV. Note that the 
largest effect comes from the shift in E^^". We have also considered a smearing in the rapidity 
of the jets, of the order of 1 unit, but this has an effect on the cross section of less than 0.1%. 
So it seems that in the longitudinal direction, the cuts of Table [1] can be directly applied at 
the parton level. 

In the following, we shall fix the splash-out correction to EJ^ ^ = O.SOEfp^ °"'^. As a 
reference we shall also choose the parametrisation given in Table [21 This is not necessarily our 
best guess, but rather one of the choices which reproduces the CDF data. We shall include 
other possibilities when we extrapolate our results to the LHC. 



6.2 Properties of the amplitude 

The accuracy of several properties and approximations presented in the literature can be 
directly tested from their effect on the cross section. 

The first property is the claim that the cross section is perturbative. We define apert as the 
value of the cross section in which all gluon momenta are larger than 1 GeV. Although this 
is not a physical observable, it is a useful quantity to test the assumptions used in theoretical 
calculations. We find that at CDF, the ratio (Jpert/<^ is 0.35 for E'^™ = 10 GeV, and falls to 
0.25 for Eip'^'^ = 35 GeV. Clearly, as already mentioned, the effect of the Sudakov form factor 
is not sufficient to shift the dominant values of the momenta into the perturbative region, and 
we see that momenta below 1 GeV contribute significantly to the cross section. Note also that 
the larger i?™", the softer the gluon loop. This is expected because at fixed proton energy, 
the production of higher-E'r jets requires larger Oj and Pi, and it is a robust property of the 
unintegrated gluon distributions that they become softer at larger x, see e.g.[M]. 

It might be, under those conditions, that modifications to the gluon propagator have to be 
taken into account [18]. It remains largely unclear however how to take these into account in 
a consistent and gauge-invariant way, as modifications to the propagator should also involve 



25 



^ 





' 




1, ^«' 1, 










10 


^•., : 


1 


i\ '■■., 


0.1 




0.01 






^'^x. 


0.001 




^--: 




nnni 







1000 



100 



> 

o 



10 



0.1 



0.01 



10 



15 20 25 
E^^'° [GeV] 

(a) 



30 35 40 



0.001 



CDF analysis 




20 40 



80 100 120 140 160 
Mjj [GeV] 



(b) 



Figure 10: (a) The effect of neglecting kj w.r.t. k or k' with the parameters of Table O (b) 
Mass distribution of the jet system, compared with CDF run II data [Ij, when changing the 
upper scale of the Sudakov form factor, as indicated. In both figures the plain curve is as in 
Fig.l 



modifications to the vertices. Also, we have not taken into account the effect of the gap 
survival factor, which will reduce the dominant values of k,. Again, the formalism to follow 
seems unclear as the protons communicate with the jet via rather soft momenta. 

The second point we can test is the neglect of the momenta kj in front of k, which we used 
in our order-of- magnitude estimate (H6l) . We show in Fig. [inia that this approximation is very 
rough, and that it overestimates the cross section at high EJ^^^. 

Finally, we confirm that the A = terms of the cross section flTSl) dominate the A = 2 
contributions, as noticed in [8]. The latter contribute 2% at ii^™" = 10 GeV to 3% at 
E^™ = 35 GeV. 

CDF also presented in [Ij the dijet mass Mjj distribution of the jet system. Note that these 
are not true data but rather predictions of ExHuMe Monte Carlo simulations [51j normalised 
to the data: the ExHuMe cross section is normalised to the data in a given Et bin, and 
the Mjj distribution comes from summing over all these Et bins. The data is presented 
for Mjj > 30 GeV, i.e. where the lower E^ cut has little influence. In Fig. (TOlb, we show 
how our calculations compare with these distributions. We have produced two curves, both 
corresponding to minimum Et of 5 GeV, and assumed that the longitudinal splash-out is the 
same as the transverse one, i.e. we took Mjj = O.SO^^sTT. The first curve corresponds to our 



reference curve. It clearly overshoots the CDF points, and predicts a considerably harder Mjj 
spectrum than the ExHuMe Monte Carlo. The second curve goes perfectly though the points. 
The only difference is that in the latter case, the upper scale of the Sudakov form factor has 
been taken as /i^ = Sgg/A, close to the choice of ExHuMe which takes /x^ = Sgg/2.62, whereas 
in the former case that scale was k2/2. We see that the two choices of the scale lead to vastly 
different result, as they affect the dependence in x = /3i//52 of the cross section. 
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Figure 11: (a) The effect of changing the upper hmit of the Sudakov loops by a factor 2, or of 
choosing Sgg/A instead of k."^; (b) The effect of changing the lower limit of the Sudakov loops 
by a factor 2. In both figures, the plain curve is as in Fig. [91 



We argued in Section 12.21 that the upper scale in the Sudakov integral for the dijet pro- 
duction should be related with the relative transverse momentum rather than the invariant 
mass of the dijet. Thus, da/dMjj obtained from experimental data without the theoretical 
bias just described would help test this point. 



6.3 Uncertainties 

So far, we have shown that it is possible, with an appropriate choice of parameters, to reproduce 
the CDF dijet data via a calculation containing several perturbative ingredients, although the 
dominant momenta are largely in the non-perturbative region. We shall now see why this is 
the case: there is no firm reason to believe the parameters of Table [21 and different reasonable 
choices can easily lead to factors of a few up or down. Basically, all the lines of Table [H can be 
changed to check on the resulting variation. We shall mention here only the most significant 
ones. We stress here that we are not trying to get the highest possible factor: our estimates 
are rather conservative, and based on changes that modest theoretical changes bring into the 
calculation. 

The main correction comes from the inclusion of vertex corrections in the form of a Sudakov 
form factor. The upper and lower limits of the integral f[T9|) can be modified while keeping 
single-log accuracy. We show in Fig. [Ill the effect of such modifications. 

As we explained above, these plots clearly show that, in this region of E-r, the vertex cor- 
rections lead to a rather uncertain situation. On the one hand, they are needed to reproduce 
the data, but on the other, their numerical impact depends on the details of their implementa- 
tion. We estimate the uncertainty coming from variations of the limits of integration as about 
a factor 3 for the lower limit, and 6 for the upper one. 
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Figure 12: (a) The differential cross section for different choices of the Sudakov form factor: 
DLA stands for the double log approximation, SEA for the splitting function approximation, 
as defined in Section 12.11 The plain curve is as in Fig. [91 (b) The uncertainty coming from 
different choices of impact factors. UgD stands for the unintegrated gluon densities, and 
LCWF for the light-cone impact factors, as explained in Section [3l 



Furthermore, we show in Fig. [T2la the effect of changes in the choice of parametrisation. 
We first show the result of a resummation of log^ terms, including angular ordering, and 
choosing the scale of as as g^. As we have seen before, the single- logs are opposite to the 
double logs, so that the cross section increases if one includes them. We also show in the same 
figure the change coming from choosing k2 as the scale of as in the Sudakov form factor. One 
sees that the choice of upper and lower limits in the Sudakov integral, and that of scale in as 
can lead to much larger effects than the inclusion of subleading logs and constant terms. 

The second uncertainty comes from the impact factor, see Fig. [12lb. Several choices of 
unintegrated gluon densities, which differ only in the choice of parametrisation of soft colour- 
singlet exchange, lead to the shaded band, while the impact factor based on a simple 3-quark 
light-cone wave function leads to the curve. Although both lead to an acceptable fit to the CDF 
data, it seems that the curvature is better reproduced by the more sophisticated unintegrated 
gluon density. Uncertainties due to possible parametrisations of the soft region, and to the 
choice of form factor, amount to a factor of at least 3. 

The parametrisations of the proton form factor used in this study contain a soft and a 
hard part. The hard part was obtained in [3ll|35] within the /cj-factorisation formalism, which 
is devised for processes with small-x gluons and not too large transverse momenta. Although 
central exclusive production takes place precisely in this regime, one can, in principle, improve 
the treatment of the hard part in the spirit of the CCFM equation [271 EH], ci-nd one can imagine 
that this improvement will shift numerical results. However, we do not expect it to reduce 
the spread of different predictions because it comes mainly from our lack of knowledge of the 
proton form factor in the soft region. 
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cuts 


A [52J 


B [53] 


«! + 02 


[0.002,0.02] 


[0.005,0.018] 


/3l+/?2 


[0.002,0.02] 


[0.004,0.014] 


\yjet 


< 1 


< 1.75 


y^luster =i|l0gfgg 


— 


<0.06 


Mn 


> 50 GeV 


> 80 GeV 



Table 3: Experimental cuts for the LHC 

To these uncertainties one must add that on the gap survival probability and that on the 
parametrisation of the splash-out, so that the overall theoretical uncertainty in the calculation 
is at least a factor 400 between the lowest and the highest estimates. 



6.4 Predictions for the LHC 

The formulae derived here can be equally used to estimate the central exclusive dijet produc- 
tion at the LHC. To do this, one needs to take into account the following corrections: 

• Specific cuts that will be used at the LHC to search for such events. These are given in 
Table [31 

• Extrapolation of the proton form factor to the LHC energies. This can be done since 
the parametrisations of the unintegrated gluon distributions are available for very small 

X. 

• Extrapolation of the gap survival probability. For an estimate, we will use the prescrip- 



tion: {SIhc) 



{s:_ 



Tevatron 



)/2. 



The predictions for the LHC based solely on theoretical calculations will be unavoidably 
plagued by the same very large uncertainty as for the Tevatron. Restricting the models with 
the Tevatron data considerably reduces this uncertainty. In Fig. [13] we show our predictions 
for the dijet central exclusive production at the LHC. Shown are two bands indicating the 
range of uncertainties. The inner band represents how different sets of parameters, tuned at 
the Tevatron to describe the central value of the data, diverge as one extrapolates from the 
Tevatron to the LHC. The outer band includes all the parametrisations presented in the other 
figures of this paper which go through all the CDF points at the la level. 

We also checked that the ratio <7pert/<7 described in Section W2\ is about 0.35 at the LHC, 
confirming that the gluon loop is still dominated by the soft region. 



7 Conclusion 

In this paper, we have seen that central exclusive production remains dominated by the non- 
perturbative region. This means that it is very important to use impact factors — such as 
the unintegrated gluon densities of [SU [351 [26] — which take the non-perturbative region into 
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Figure 13: The reference curve extrapolated to the LHC, for the two sets of cuts shown in 
Table [3l The graph to the left corresponds to the cuts given in column A and the one to the 
right to cuts B. The inner bands give the theoretical error for curves reproducing the central 
values of the CDF data, and the outer bands correspond to the 1 a errors. 



account. We have shown that this dominance of the soft region implies large uncertainties in all 
the ingredients of the calculation: vertex corrections, impact factors and screening corrections. 
We evaluate the uncertainty as being a factor 20 up or down, with no theoretically preferred 
curve. 

At present, one can tune a perturbative calculation to the CDF run II data on dijet 
exclusive production, and try to use it to predict the cross sections for the production of other 
systems of particles. 

There are several problems with this approach. The first one is that the huge Sudakov 
form factors of the dijet case suppress diagrams where the hard scale is concentrated in one 
propagator of the graph. Other graphs, such as those of Fig. [TJ where the hard scale flows 
through the diagram, are suppressed by propagators, but one may expect the vertex corrections 
to be substantially smaller. So using dijet production as a handle on e.g. Higgs production 
may be misleading. 

The second problem concerns the extrapolation to LHC energy. This is due to the fact that 
unitarisation effects will be important when going from 2 TeV to 14 TeV. The embedding of 
the process in a multi-pomeron environment can be achieved for short-distance hard partonic 
processes, leading to the usual gap survival probability. However in general, if the partons are 
not at very small distances, the gap survival probability will be more complicated, and larger. 
Also, it is usually assumed that Regge factorisation holds, i.e. that the process can be written 
as a product of factors depending on ti = fc^ and t2 = k^, which is not the case here. 

Clearly, all the above questions can only be settled by comparing to more data. Hence we 
believe that a measurement of quasielatic jet production at the LHC, including the rapidity 
distribution, the mass distribution, as well as the Et distribution, will be a very important 
task that will help to clarify many of the issues raised in this paper. 
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Appendices 

A Lowest-order two-gluon production 
in the quasi-multi-Regge kinematics 

Consider the production of a colour- singlet two-gluon state in the qq collision. The two gluons 
are not required to have large transverse momenta or large invariant mass, but they are well 
separated in rapidity from the quarks. Besides, we require as usual that there is no overall 
colour flow in the t-channel. Then, the kinematics and all the simplifications discussed in 
Sect. 11.11 and 11.21 still hold, and we are left with two generic diagrams shown in Fig. |H 

From the BFKL point of view, we deal with two-gluon production in the quasi-multi-Regge 
kinematics (qMRK). The imaginary part of the amplitude can be written as 

„6 Xab r 

ImM = ImMa + lmMb = s^-— d^ke'i^^e'-'^^* (57) 



„6 xab r 


-, 


_k2(k + ki)2(k + k3)2 k2(k + ki)2(k + k2)2(k + k^ - 


-kapj 



Here C^^ and C^^ are the usual Lipatov effective vertices (the tilde indicates that this vertex 
originates from the other t-channel leg), e.g. 

C^^ = P"' (a, - 2^^±^] - q^^ (p, - 2^^±^] + (2k + k^ + k^)'^^ , (58) 

while €2^^'^ is the effective RRGG vertex for two gluon production in qMRK [22]. In the multi- 
Regge kinematics (MRK) limit, /3i ^ /32, it factorizes into usual Lipatov vertices describing 
successive production of two MRK gluons: 



On the other hand, when |k2| ^ |ki|, |k3|, one would recover from (1571) the amplitude ([8]) 
described in the main text. 

The two terms in flFTI) describe the respective contributions of the two generic diagrams of 
Fig. m Let us estimate how the contribution of the second diagram compares with the first 
one in the case of large \i^ and, for simplicity, at [3i ^ /32- 
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The square of the first diagram is 
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For nearly forward scattering, |ki|, \k^\ <^ |k|, |k'|, so that the estimate simplifies to 
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A similar estimate for the second diagram gives 
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which for the forward case becomes 
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Thus, the second diagram is not only suppressed with respect to (I6T]) by an extra power of 
k|, but also enhanced by an extra logarithm. If the loop corrections to this diagram can be 
resummed and presented in a form analogous to the Sudakov form factor, then the k2-integral 
will be shifted towards large values and no logarithmic enhancement in this new suppressing 
factor will occur. Hence the new Sudakov form factor may be much larger than the usual one. 

B Intermediate production of gluons in multi-Regge 
kinematics as a source of corrections 
to central exclusive production 

If a colourless high-mass system / can be "radiated" off the t-channel gluons, then the same 
system can be produced via two intermediate gluons in multi-Regge kinematics (MRK) fol- 
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Figure 14: a: Generic diagram representing production of two gluons in MRK followed by their 
fusion into the high mass system, h: extra gluon rungs above or below the central subprocess 
can be incorporated into the definition of the unintegrated gluon density, c: an extra gluon 
rung that couples the central subprocess with the screening t-channel gluon represents a single- 
log correction to the vertex, d: higher-order single-log enhanced correction to diagram Fig.[Hlc. 



lowed by their fusion into /, see Fig. El We do not present here a detailed calculation of this 
process; however, it can be immediately seen from (jl]) and ( 157]) that the resulting cross section 
is flat in rapidities of the intermediate gluons. Integrating over the corresponding lightcone 
variables but keeping M? fixed yields an extra longitudinal logarithm. The particular dia- 
gram shown in Fig. [THa is in fact a double-log-enhanced correction to the direct production 
of system /. In the standard calculation it is effectively included via the Sudakov form factor. 

Even higher-order corrections to this diagram include diagrams with extra gluon rungs that 
join the two t-channel gluons. If these extra gluons are above or below the central subprocess, 
as in Fig. (TUb, they can be absorbed into the definition of the unintegrated gluon density. 
However, if the gluon is attached inside the central subprocess. Fig. [Hlc, it leads to a correction 
to the effective vertex and cannot be absorbed into the gluon density. 

This correction is enhanced by a single logarithm. If the lightcone variables of the two 
intermediate gluons are /3i ^ (32, then the logarithm that appears is of type 



log 



P2 



log 



Ml 
q2 



where q^ is the (moderate) transverse momentum of the intermediate gluons. Thus, this dia- 
gram represents a single-log enhanced correction to the central exclusive production of system 
/, which involves the screening gluons and which is absent in the pure gg ^ f process. Since 
the single-log enhanced terms in the Sudakov integral are very important for the kinematic 
range considered, one must take seriously such corrections. 

As the production of two gluons in MRK requires gluon emission from both t-channel legs, 
one encounters additional single-log enhanced diagrams such as shown in Fig. [TUd. 

We would like also to stress that, even without the above estimates, there is a need to 
evaluate such "non-conventional" diagrams, which stems directly from self-consistency of one's 
approach. Indeed, if one believes the diagrams such as Fig. (TUd are not logarithmically 
enhanced, then the bare diagram Fig. [Hlc does not get suppressed at all. It can then become 
comparable with the usual Sudakov-suppressed diagrams, as was explained in Section 12.31 If, 
on the other hand. Fig. [TUd and similar diagrams are logarithmically enhanced, one needs to 
resum them to estimate their effect. 
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